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Abstract. In an effort to extend the still limited sample of s-process enriched post-AGB stars, high-resolution, high signal- 
to-noise VLT+UVES spectra of the optical counterparts of the infrared sources IRAS 06530-0213 and IRAS 08143-4406 were 
analysed. The objects are moderately metal deficient by [Fe/H] = -0.5 and -0.4 respectively, carbon-rich and, above all, heavily 
s-process enhanced with a [ls/Fe] of 1.8 and 1.5 respectively. Especially the spectrum of IRAS 06530-0213 is dominated by 
transitions of s-process species, and therefore resembling the spectrum of IRAS 05341+0852, the most s-process enriched object 
known so far. The two objects are chemically very similar to the 21pm objects discussed in Van Winckel & Reyniers (2000). 
A homogeneous comparison with the results of these objects reveals that the relation between the third dredge-up efficiency 
and the neutron nucleosynthesis efficiency found for the 21yum objects, is further strengthened. On the other hand, a detailed 
comparison with the predictions of the latest AGB models indicates that the observed spread in nucleosynthesis efficiency is 
certainly intrinsic, and proves that different 13 C pockets are needed for stars with comparable mass and metallicity to explain 
their abundances. 

Key words. Stars: AGB and post-AGB - Stars: abundances - Stars: carbon - Stars: individual: IRAS 06530-0213 - Stars: 
individual: IRAS 08143-4406 



' 1 . Introduction There are only few sources displaying a spherically symmetric 

morphology, most sources possessing a clear bipolar structure 

£j ■ In the last two decennia, the progress in the theoretical mod- (see e g Sahai 2001) The aspherical structure was already 

Ctf ; elling of AGB stars, as well as the qualitative and quantitative demori strated earlier by ground based observations both in op- 

tmprovement of the observational data of (post-)AGB stars, are tical images of me scattered visible light (e _ g _ Hrivnak et al. 

impressive. Whereas, before the launch of the IRAS satellite, l99g ^ R ^ ."75), and in mid-IR images of the thermal emis- 

only a few post-AGB candidates were known, the present sam- sion of dust (e g Meixner et al . 1999? R „ V >y Many more 

pie consists of about 220 objects (Szczerba et al. 2001). The stmctural details werej however, resolved on HST images (e.g. 

main keyword, however, that can be applied to the post-AGB Ueta et al 2000; Hrivnak et aL 20 01). For a recent review, we 

sample is diversity. In practically all their aspects, post-AGB re f er j G Balick & Frank (2002) 
stars are much more diverse than theoretically anticipated. 

T , , ii . . . Concerning the circumstellar chemistry, the high resolution 
In the case of the mo rphology, the resolved post-AGB stars b J b 
,., - i j »fu spectroscopy of the infrared excess also reveals an interest- 
display a surprisingly wide variety of shapes and structures. r rj 

ing chemical evolution. The observations made by the Infrared 

Space Observatory (ISO) satellite brought important progress 
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„„.„,-, , in this domain. In an O-nch chemistry, the dust is mainly com- 
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* t, j u ii t a t <u c c >u posed of amorphous silicates, but also bands of crystalline sil- 

Based on observations collected at the European Southern r r > j 

Observatory, Paranal, Chile (ESO Programme 66.D-0171) icates 818 found - In the carbon dominated chemistry, the dust 

" Table 4 is only available in electronic form at the CDS via consists mainly of amorphous carbon, SiC, and probably MgS 

anonymous ftp to cdsarc.u-strasbg.fr (138.79.128.5) or via (the carrier of the broad 30yum band, Hony et al. 2002). Also the 

http : //cdsweb . u-strasbg . fr/cgi-bin/qcat?l/A+A/ "Unidentified InfraRed" (UIR) bands are seen in these IR spec- 
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(PAHs). A sub-class of C-rich post-AGB stars shows a feature 
at 21pm, discovered by Kwok et al. (1989), possibly attributed 
to TiC nanocrystals (von Helden et al. 2000). 

Finally, there is an intriguing diversity in the observed pho- 
tospheric chemical patterns of the post-AGB stellar sample. 
The third dredge-up is expected to bring carbon and freshly 
synthesized s-process elements to the photosphere of an AGB 
star. The spectral lines of these elements are thus expected to 
be strong in the spectra of post-AGB stars. This is indeed what 
we observe, but remarkably only in a certain sub-class of post- 
AGB objects. There are very similar objects (with a compa- 
rable mass and metallicity) that do not show these enhance- 
ments. This chemical dichotomy is surprisingly strict in the 
sense that an object is either severly enriched, or it is not en- 
riched at all (even s-process deficiencies are observed). This is 
illustrated in Fig. 3 of Van Winckel (2003) on which the over- 
abundance of the s-process element zirconium (atomic number 
Z = 40) is plotted against metallicity for a sample of stars that 
all show similar characteristics for a genuine post-AGB classi- 
fication, like a double peaked SED, low metallicity and kine- 
matics pointing to membership of an old population. The chem- 
ical difference is directly clear from the spectrum itself and we 
refer to Fig. 1 for an illustration. Moreover, all s-enriched ob- 
jects studied till now are suspected or confirmed 21pm sources. 
Six of these sources were recently analysed in a homogeneous 
abundance study by Van Winckel & Reyniers (2000), hereafter 
Paper I. In the non s-enriched objects (see e.g. Luck et al. 1990; 
Van Winckel 1997, for abundance analyses of these objects), 
no enrichments were found larger than the [Zr/Fe] ~ 0.2 error 
level. 

In the present study, we present the analysis of two newly 
discovered s-process enriched post-AGB stars, IRAS 06530- 
0213 and IRAS 08143-4406 (Tab. 1). Despite the fact that 
these IRAS sources have been both candidate post-AGB stars 
for more than a decade (as derived from their position in the 
IRAS color-color diagram, e.g. Preite-Martinez 1988), they are 
both still poorly studied. For IRAS 06530-0213 this can partly 
be attributed to the optical weakness of this source (Tab. 1), 
while for IRAS 08143-4406 it is less clear why this source has 
been neglected for such a long time. Based on a low resolu- 
tion spectrum (10.7 ApixeL 1 ) taken with the 2.5 m INT on La 
Palma, Slijkhuis (1992) attributed a spectral type F0I to the op- 
tical counterpart of IRAS 06530-0213. In a very recent study, 
Hrivnak & Reddy (2003) found that this is a too early spec- 
tral type and classified IRAS 06530-0213 as a F5 supergiant. 
The spectral type of IRAS 08143-4406 (F8I) was taken from 
Reddy & Parthasarathy (1996) and is based on a low resolution 
(5.7 A pixel -1 ) spectrum taken with aim telescope. 

The paper is organised as follows: In Sect. 2 we describe the 
observations of the two programme stars, together with the data 
reduction. In Sect. 3 the radial velocities of the two programme 
stars are compared with literature values, with the aim of test- 
ing their possible binary nature. Sect. 4 focusses on the tech- 
nical aspects of the analysis, including atmospheric parameter 
determination, line selection and spectrum synthesis. Also our 
abundance results are presented in this section. In the next sec- 
tion, we compare the two program stars with the 21yum stars of 
Paper I. It turns out that the two stars share the same chemical 



Table 2. Observational log. Spectral gaps occur between 577 nm and 
583 nm and between 854.4 nm and 864.5 nm due to the spatial gap 
between the two UVES CCDs. 



date 


UT 


exp.time 


wavelength 


S/N 




start 


(sec) 


interval (nm) 




IRAS 06530-0213 


2001-01-12 


03:36 


3x1800 


477.5-681 


100 


IRAS 08143-4406 


2001-01-16 


06:41 


2x1800 


374.5-498 


130 


2001-01-16 


06:41 


1x1800 


670.5-1055 


sat. 


2001-01-16 


07:16 


3x500 


670.5-1055 


200 


2001-02-01 


04:21 


1x1800 


477.5-681 


150 



sat. = saturated 

properties with the 21pm stars. Sect. 6 is devoted to a detailed 
comparison with theoretical chemical AGB model predictions. 
In this section, also the 21pm stars are compared with these 
chemical models. Finally, in Sect. 7 we summarize our most 
important findings. 

2. Observations 

High resolution, high signal-to-noise VLT+UVES spectra of 
the two programme stars are taken in the framework of our on- 
going program to study the photospheric chemical composition 
of stars in their last stages of evolution (e.g. Paper I; Reyniers 
& Van Winckel 2001, 2003). The two targets are part of a larger 
sample of eleven post-AGB objects that were observed in ser- 
vice mode during two periods (04-07/2000 and 01-02/2001, see 
Reyniers 2002). The resolving power of these spectra varies be- 
tween ~55,000 and ~60,000. Some details about the observa- 
tions are given in Tab. 2. 

The standard reduction was performed in the dedicated 
"UVES context" of the MIDAS environment and included bias 
correction, cosmic hit correction, flat-fielding, background cor- 
rection and sky correction. We used optimal extraction to con- 
vert frames from pixel-pixel to pixel-order space. The spec- 
tra were normalised by dividing the individual orders by a 
smoothed spline function defined through interactively identi- 
fied continuum points. Orders were merged after this normali- 
sation. For a detailed description of the reduction procedure, we 
refer to Reyniers (2002). In Tab. 2, we also list some indicative 
signal-to-noise values of the final data product. Sample spectra 
can be found in Figs. 1, 3 and 5. 

3. Radial velocities 

Since binarity can influence the chemical evolution of a post- 
AGB star drastically (see e.g. Van Winckel et al. 1995), it is 
important to determine and monitor accurate radial velocities 
for these objects. The large wavelength coverage together with 
the high resolution of our UVES spectra permit us to deter- 
mine these velocities using a large number of lines. The ob- 
tained radial velocities are gathered in Tab. 3, together with the 
few values found in the literature. An additional velocity is ob- 
tained with the CORALIE spectrograph mounted on the Swiss 
telescope Euler in La Silla, Chile. Highly accurate velocities 
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Table 1. Basic parameters of the two objects discussed in this study. 
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IRAS 


Equatorial 


Galactic 


Visual 


Spectral 




IRAS Fluxes (Jy) 






coordinates 


coordinates 


magnitude 


Type 


fn 


fa 


feo 


/too 




Q '2000 


^2000 


1 


b 


m(b) 


m(v) 












06530-0213 


06 55 32.1 


-02 17 30 


215.44 


-0.13 


16.4" 


14.0" 


F0I" 


6.11 


21 Al 


15.05 


4.10 












16.3* 


14.1* 


F5F 










08143-4406 


08 16 02.9 


-44 16 01 


260.83 


-5.07 


14.1* 


12.4* 


F8I* 


0.60 


9.26 


6.06 


<3.73 



Source SIMBAD, except: " Slijkhuis 1992, * Reddy & Parthasarathy 1996, c Hrivnak & Reddy 2003 




zoo 

Is 

■-- hs 



312 612.5 613 

Wavelength (nm) 



Fig. 1. The spectra of IRAS 06530-0213 and IRAS 08143-4406 compared with the heavily s-enriched post-AGB star IRAS 05341+0852 and 
the non s-enriched post-AGB star HD 108015. IRAS 05341+0852 is discussed in Paper I, and it is the most heavily s-process enriched star 
known so far. HD 108015 has been analysed by Van Winckel (1997) and has atmospheric parameters comparable with the other three (T e g, 
logg, £,) = (7000 K, 1.5, 4.0km s~') but is slightly less metal deficient ([Fe/H] = -0.1). A complete line identification of this spectral interval 
has been pursued by the use of the VALD database (Kupka et al. 1999). Lines of light s-process elements (Sr peak) are identified by a dotted 
line; lines of heavy s-process elements (Ba peak) by a dash-dotted line; lines of other elements (which are mainly a and iron peak elements) 
by a full line. The s-process enhancement of the three IRAS stars is clear. IRAS 05341+0852 is the most s-enriched one, but the stronger lines 
are also due to a slightly lower temperature of this object. All spectra are VLT+UVES spectra, except the spectrum of HD 108015, which is a 
ESQ1.5+FEROS spectrum taken at March 22, 2000. 



are obtained with this instrument by cross-correlation observed 
spectra with a spectrum mask of comparable spectral type. 

It is obvious that there is still monitoring needed to draw 
any final conclusion, but till now, there is no evidence for a 
binary motion for the two programme stars. 



4. Abundance analysis 

4.1. Atomic data and atmospheric parameters 

A list of lines useful for the chemical analysis of A and F type 
stars has been collected at the Instituut voor Sterrenkunde dur- 
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Table 3. Heliocentric radial velocities of the programme stars, both 
from the literature and from this study. If the velocity is obtained in 
this study, the number of lines on which it is based is given in the last 
column; otherwise the reference is given in the same column. The last 
velocity is obtained with the CORALIE spectrograph. Typical errors 
on these velocities are between 1.0 and 1.5 km s -1 . 



Date 




method 1 


ref.* or number 


(yyyy-mm-dd) 


(kms 4 ) 




of used lines 


IRAS 06530-0213 


1991-04 


50 


a 


ref. 1 


1997-10-17 


51.0 


b 


ref. 2 


2001-01-12 


49.8 


b 


n = 236 


2001-12-10 


50.4 


b 


ref. 2 


IRAS 08143-4406 


2001-01-16 


51.5 


b 


n= 134 


2001-02-01 


52.2 


b 


n = 434 


2002-05-22 


49.6 


c 





''method a: CO (J = 2-1) line; b: mean of individual optical lines 
c: cross-correlation with spectrum mask t ref. 1: Hu et al. 1994; 2: 
Hrivnak & Reddy 2003 



Table 4. Line identifications and measured equivalent widths for 
IRAS 06530-0213 and IRAS 08143-4406. This table is only available 
in electronic form at the CDS. It contains the following columns: the 
rest wavelength (in A), the identification, the lower excitation poten- 
tial (in eV), the adopted log gf value, the measured equivalent widths 
for IRAS 06530-0213 (in mA), the measured equivalent widths for 
the red setting (2001-02-01, see Tab. 2) of IRAS 08143-4406 and the 
measured equivalent widths for the dichroic setting (2001-01-16, see 
Tab. 2) of IRAS 08143-4406. 



ing the past few years and is regularly updated. This list is de- 
scribed in detail in Paper I. A recent update of the list after the 
publication of Paper I consists of the inclusion of the data pub- 
lished by Lawler et al. (2001a) for La and Lawler et al. (2001b) 
for Eu. Using this same list, we do not only restrict ourselves 
to lines with reliable atomic data, but we also ensure that the 
different analyses are perfectly homogeneous. 

Equivalent widths of unblended lines were measured with 
direct integration; multiple gaussian fitting was applied for 
lines with (a) blended wing(s). The entire line list, together with 
the measured equivalent widths for both stars can be found in 
Tab. 4, which is published in electronic form at the CDS. The 
atmospheric parameters were obtained by the usual spectro- 
scopic method in which the effective temperature is obtained 
by imposing the iron abundance, derived from the individual 
Fe i lines, to be independent of lower excitation potential; the 
gravity by imposing ionization equilibrium; the microturbulent 
velocity by imposing iron abundance derived from the individ- 
ual Fe i lines to be independent of (reduced) equivalent width. 
The model atmospheres of Kurucz (1993) were used, in com- 
bination with the latest version (April 2002) of the LTE abun- 
dance calculation routine MOOG (Sneden 1973). 

The results on other species can be used as a check for the 
atmospheric parameters, providing that the number of lines for 
that species is large enough. Elements for which we detected 



lines of different ions can be used as a first consistency check 
of the model atmosphere. A temperature check is difficult since 
for most elements, the range in excitation potential is quite 
small. Therefore, only the microturbulent velocity parameter 
ft can be checked in this way. For a given element, the most 
straightforward method is then to search for that f t for which 
the slope of the linear fit on a reduced width - abundance di- 
agram equals zero. With this method, however, one loses the 
particular characteristics of the specific diagram, in the sense 
that mostly outliers in log(W^//l) determine the obtained slope. 
Therefore, it is often much more instructive to inspect the di- 
agram by eye, than to rely on a least-squares method. This is 
done in Fig. 2, on which we plotted the Fe, Ce and Nd dia- 
grams for IRAS 06530-0213, and the Fe, La and Nd diagrams 
for IRAS 08143-4406. From this figure it is clear that (a) the 
spread in reduced width log(W / i//l) is significantly larger for Fe 
than for the other elements and that (b) the spread in abundance 
log e is smaller in the case of iron. For the non-iron species on 
Fig. 2 it is difficult to decide whether a real trend is seen or not. 
Especially the points in the Nd diagram display a somewhat 
confusing pattern. Several authors already stressed the need for 
new atomic data for Nd n (see e.g. Sneden et al. 2002, for a dis- 
cussion of the available Nd n data). We can conclude that there 
is no clear indication for a dependence of f t on the species that 
is used to derive it. Note that an increase of ~3 km s _1 does not 
induce large changes in the derived abundances (a decrease of 
~0.1 dex for the strongest lines in our line list). 

A particular problem arose in the determination of the pa- 
rameters of IRAS 08143-4406. Contrary to IRAS 06530-0213, 
the star is measured in two different spectrograph settings: on 
2001-01-16 the dichroic setting (hereafter DI2) is used with 
a simultaneous observation in the blue and in the red, while 
on 2001-02-01 the red setting centered on 580nm (hereafter 
RED580) is observed (see also Tab. 2). Between these two ob- 
servations, there is a time gap of 16 days. By inspecting the 
spectral overlap of the two observations, we concluded that the 
atmospheric parameters apparently changed during this time 
gap. The adoption of two different sets of atmospheric param- 
eters for each setting is therefore necessary. Unfortunately, the 
number of Fe lines suitable for parameter determination is too 
low for the dichroic setting. An alternative method had to be 
developed to obtain the parameters for this setting. 

We chose to deduce the atmospheric parameters for the 
2001-01-16 setting of IRAS 08143-4406 relative to its param- 
eters of the 2001-02-01 setting, by demanding that the lines in 
the overlap should yield the same abundance. So, after determi- 
nation of the atmospheric parameters for IRAS 08 143-4406 for 
the 2001-02-01 setting, we searched for the model parameters 
for the 2001-01-16 setting, by minimizing the abundance dif- 
ference of the lines in the overlap. This yielded model parame- 
ters (T eS , logg, = (7250 K, 1.5, 5.5 kms" 1 ) for the 2001-02- 
01 setting, and (7050 K, 1.2, S^kms" 1 ) for 2001-01-16. One 
should note that a difference of 0.3 dex in gravity is rather large: 
assuming a mass of 1 .5 M Q , a log g of 1 .5 yields an approxima- 
tive radius of R ^ 36 R Q , while a log g of 1 .2 yields R ^ 5 1 R . 
This results in a pulsation velocity of v^7.5kms _1 , which is 
not consistent with the seemingly constant radial velocity of 
IRAS 08143-4406 (see Tab. 3). Moreover, such a change in at- 
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Fig. 2. Reduced equivalent width - abundance diagrams for elements 
of which the abundance is based on more than fifteen lines. Such dia- 
grams are constructed to derive the microturbulent velocity for the 
model atmosphere; the correct value for £ t is found if no trend is seen 
on the diagram. In this figure, we check if the abundances of some 
non-iron species are consistent with the value for f, that was derived 
from a study of the iron lines. 

mospheric parameters would induce a difference of ~0.6 mag 
in V in 16 days. Since the gravity is not well constrained (a 
typical error in log g is 0.5), we decided to keep the difference 
in the parameters in order to keep the abundances deduced from 
the two settings as consistent as possible. 

4.2. Spectrum synthesis of IRAS 06530-021 3 

4.2.1 . The macroturbulent broadening g m 

When computing a spectrum synthesis, a broadening parame- 
ter, additional to the usual atmospheric parameters, is required 
to match the synthetic line profiles with the observed ones. This 
broadening is the combined effect of instrumental, macrotur- 
bulent and rotational broadening. The latter one is, however, 
thought to be small in F-G supergiants. Since the instrumental 
broadening can be deduced from the width of the Th-Ar lines 
of the calibration spectra, the macroturbulent broadening (£ m ) 
is the only parameter that has to be estimated before comparing 
synthetic spectra with observed ones. We concentrated on the 
spectra of IRAS 06530-0213, as this star is the most interesting 
in the framework of s-process nucleosynthesis. 
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Fig. 3. Synthesis of two lines of neutral calcium. From this figure, it is 
clear that spectral lines cannot be fitted using one fixed value for the 
macroturbulent broadening f m . The synthesis is made using the abun- 
dance as derived from their equivalent width. The two abundances are 
in good agreement. The best-fit macroturbulent broadening, however, 
differs by 5.5kms~'. 

For the instrumental broadening we took the median of the 
resolution of the Th-Ar lines used in the calibration, which is 
5A = /1/59525 A for the lower (EEV) CCD and 6A = i/56575 A 
for the upper (MIT) CCD. The macroturbulent velocity % m 
was then obtained as follows. We selected 121 unblended lines 
from the abundance analysis of IRAS 06530-0213 from differ- 
ent species and ionisation (neutral as well as singly ionised). 
Each of these 121 lines was synthesized with MOOG using the 
abundance as derived by its equivalent width. The macrotur- 
bulent broadening % m is the only free parameter in this synthe- 
sis. We fitted the synthetic spectra by varying this parameter 
in steps of 0.5kms~'. The profile used for the macroturbulent 
broadening is a radial-tangential macroturbulence profile based 
on the work of Gray (1992). During this fitting procedure, we 
soon realized that we could not fit the lines with one fixed value 
for % m , but that each line required its own value. This is illus- 
trated in Fig. 3, on which we synthesized two calcium lines. 
These lines are lying in the same spectral region and the abun- 
dances derived from their equivalent widths are in very good 
agreement. However, we need two different values for £ m in a 
synthesis of their profiles. 

Analyzing these different values for t; m , we found a sur- 
prisingly tight correlation of the macroturbulent velocity with 
the (reduced) equivalent width (classical correlation coefficient 
p = 0.81). A simple least squares fit gives 

£ m = 10.9 log(W,/i) + 71.7 

Further inspection of the results gives an offset of ~2.3kms~ 1 
between lines of neutral and lines of ionised species with the 
same reduced width. Therefore, we can specify the previous 
relation by making the distinction between lines of neutral and 
lines of ionised species: 

^(neutral lines) = 10.9 log(W^M) + 70.2 ± 2.7 
^(ionised lines) = 10.9 log(W^M) + 72.4 ± 3.2 
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Fig. 4. Relation between the (reduced) equivalent width logCWi/zi) 
and the macroturbulent velocity £,„ of 121 lines used in the analysis 
of IRAS 06530-0213. The lines from ionised species have systemati- 
cally higher macroturbulent velocities than lines from neutral species 
of the same width (upper panel). Therefore, two separate relations are 
deduced (middle and lower panel). The errorbars are the standard de- 
viations on the mean trend multiplied by 1.5. 



where the errorbar is empirically deduced from Fig. 4 (being 
1.5 times the standard deviation on the mean trend). The cor- 
relation coefficient for lines of neutral species is p n = 0.90; the 
same coefficient for lines of singly ionised species is p; = 0.70. 
The reason for this dependence is not clear, nor a description 
of this effect was found in literature. This might be due to 
the fact that most authors deduce the macroturbulent velocity 
from just a few lines or lines in the immediate proximity of 
the line/region under interest. We interpret this effect as an op- 
tical depth effect of the macroturbulent velocity, the stronger 
lines being formed on average at lower geometrical depth. We 
want to stress that a similar relation, however less documented, 
was found for HD 172481 (Reyniers & Van Winckel 2001), 
so it is certainly not an isolated peculiarity of the spectrum of 
IRAS 06530-0213 studied here. 

4.2.2. O synthesis 

Due to the high excitation potential of the oxygen lines, a re- 
liable oxygen abundance is always difficult to derive in this 
temperature-gravity domain, and mostly only the oxygen triplet 
at 6155A is accessible. The situation is even worse for heav- 
ily s-enriched objects like IRAS 06530-0213 since the triplet is 
significantly blended by lines of s-process species. Important 
blends at the position of the triplet are lines of Fei (6157.73 A), 
Prn (6157.81 A) and Ndn (6157.82 A). Therefore, the only 
way to obtain a reliable oxygen abundance in heavily enriched 



Table 5. Results of the spectrum synthesis of the Oi triplet of 
IRAS 06530-0213. 



A 


X 


logg/ 


adopted f m 


1°£ 




(A) 


(eV) 




(kms- 1 ) 






6155.971 


10.74 


-0.674 


13.5 (±3) 


8.58 (d 


=0.06) 


6156.778 


10.74 


-0.453 


15.25 (±3) 


8.71 (d 


=0.09) 


6158.187 


10.74 


-0.307 


16.25 (±3) 


8.62 (d 


=0.05) 










8.64 





objects is by spectrum synthesis. For the O lines themselves, 
we used the atomic data published by Biemont et al. (1991). 

Furthermore, the relation between equivalent width and 
macroturbulent velocity presented above causes a double diffi- 
culty when making the synthesis. First, the relation is deduced 
for single, non-blended lines. Hence, it is not clear how the 
relation should be applied to blended lines. Second, since the 
three lines are of different strength, each line requires its own 
macroturbulent broadening, so that the lines of the triplet had 
to be separately fitted. As a consequence, it is also difficult to 
show the synthesis on one single figure. Therefore, the results 
are summarized in a table (Tab. 5). 

4.3. Elements beyond the Ba peak 

Inspired by the detection of a hafnium (Hf, Z=72) line in the 
objects IRAS 07134+1005 and IRAS 19500-1709 as reported 
in Paper I, we initiated a systematic search for other elements 
beyond the Ba peak. Our VLT+UVES spectra of the heavily 
enriched objects are ideally suited for this search since they 
combine a large wavelength coverage with a high signal-to- 
noise. Moreover, the release of a new database with atomic 
data of the lanthanides (Z = 57-71) at Mons University (B), 
the so called D.R.E.A.M. project, is an invaluable supplement 
to the VALD database which was our only source of atomic 
data of species in this mass range till now. The potential of 
this new database was recently illustrated by Reyniers et al. 
(2002). Since it also contains atomic data for doubly ionised 
lanthanides, we extended our search to these ions. The details 
of our successful search are published in a dedicated letter 
(Reyniers & Van Winckel 2003), while here we only report 
on the results of the two programme stars. We have derived 
abundances of gadolinium (Gd, Z=64), ytterbium (Yb, Z=70), 
lutetium (Lu, Z=71) and tungsten (W, Z=74) for IRAS 06530- 
0213, and of Gd and Lu for IRAS 08143-4406. 

4.4. Hyperfine structure 

In Paper I we studied the influence of hyperfine splitting (hfs) 
on the abundances of elements which are considered to be sen- 
sitive for this effect. The general conclusion of this exercise 
was that the hfs effect on the lines used in the abundance anal- 
ysis is very small, due to the relatively high effective tempera- 
ture of the 21yum stars. This is only true for weak lines, while 
for stronger lines an effect is indeed expected. Such strong 
lines were, however, not used in the abundance determination. 
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Fig. 5. Synthesis of the 6221.890 A Lu line in IRAS 06530-0213, with- 
out (left panel) and with (right panel) hyperfine decomposition. The 
adopted abundance is indicated below. The difference in abundance 
illustrates the necessity of the inclusion of hfs for this line. 

Because several papers discussing hfs constants for different 
elements recently appeared in literature, we decided to repeat 
this exercise for the two program stars in this paper. The influ- 
ence of hfs was determined by calculating the equivalent width 
of a studied line, with and without hfs decomposition. In this 
way, the exact profile of the line is eliminated, and no broad- 
ening factors have to be applied. We focussed on the s-process 
elements La, Eu and Lu. 

For La, hfs constants A and B were taken from Lawler 
et al. (2001a). The La line list of IRAS 06530-0213 contains 
only one line with hfs constants for both levels (5482.268 A, 
W A =91.8 mA). The effect of hfs on the abundance is smaller 
than 0.01 dex. The La line list of IRAS 08143-4406 contains 
eight lines for which there are hfs constants in Lawler et al. 
(2001a). For seven of these lines, the effect is negligible 
(<0.02 dex), while for the strongest line in the list (6262.287 A, 
W A = 148.0mA), the effect is significant. An hfs-treatment of 
this line yields an abundance of 2.60 dex, which is 0.15 dex 
lower than a treatment without hfs. 

The Eu abundance is often considered to be strongly hfs 
sensitive, and an hfs treatment is indispensable when study- 
ing Eu abundances in cool stars. We studied the strongest Eu 
line of the two stars that was used in our analysis (6645.064 A, 
103.2 mA in IRAS 06530-0213). A and B constants for this line 
were taken from Lawler et al. (2001b). The solar isotopic com- 
position was applied. It turned out that, even for this relatively 
strong line with a strong decomposition, the effect on the de- 
rived abundance is marginal (<0.02dex). We also studied the 
effect of hfs on the Eu abundance of the 21//m stars of Paper I. 
Only for the cooler and strongly enriched IRAS 05341+0852, 
the effect was noticeable (an abundance decrease of 0.05 dex 
and 0.07 dex for the lines at 6437.640 A and 7194.830 A re- 
spectively). 

While for La and Eu, the effect of hfs on the abundance 
is much smaller than other uncertainties (loggf, continuum 
placement, undetected blends, . . . ), the situation is different 



for the Lu abundance. Hfs constants for the two detected Lu 
lines were taken from Brix & Kopfermann (1952) and Den 
Hartog et al. (1998). In the case of IRAS 06530-0213, the dif- 
ference in abundance between a hfs and a non-hfs synthesis is 
-0.27 dex for the 6221. 890A line (see Fig. 5), and -0.14dex 
for the 6463.107 A line. For the less s-enhanced IRAS 08143- 
4406, hfs only affects the line profile of the two Lu lines, but 
not their equivalent width. 

4.5. Abundance results 

In Tab. 6 we present the complete abundance analysis results 
of IRAS 06530-0213 and IRAS 08143-4406. The same results 
are graphically presented in Fig. 6. On this figure, the different 
groups of elements are marked with different symbols. We will 
now summarize the main results for each of these groups. 

Metallicity Both stars are (moderately) metal deficient, 
with iron abundances of [Fe/H] = -0.5 and -0.4 for 
IRAS 06530-0213 and IRAS 08143-4406 respectively. The 
other iron peak elements all follow this deficiency, al- 
though zinc (Zn) is somewhat higher than expected for 
IRAS 08143-4406. 

CNO-elements Both stars are clearly carbon enriched, with 
a huge enhancement in the case of IRAS 06530-0213 
([C/Fe] = +1.0). As a consequence, we derive a high C/O 
number ratio of C/O = 2.8 for this star; For IRAS 08143- 
4406, this ratio is 1.3. One has to note, however, that the 
errors on mainly the oxygen abundance preclude an accu- 
rate C/O number ratio determination based on photospheric 
atomic lines. 

a-elements The simple mean of the [el/Fe] values of the 
(available) a-elements Mg, Si, S, Ca and Ti yields for both 
stars [a/Fe] = +0.2. Such an enhancement is normal for 
stars in this metallicity range, as a consequence of the 
galactic chemical evolution and therefore does not corre- 
spond to an intrinsic enhancement. 

s-process elements From Fig. 6, it is clear that the s-process 
enrichment of the two objects under study is very strong. 
This enrichment is the most important argument for the post 
third dredge-up status of the two stars. 

The s-process elements observed in evolved stars can be di- 
vided into two groups: the light s-process elements around the 
magic neutron number 50 (Sr, Y, Zr) and the heavy s-process 
elements around the magic neutron number 82 (Ba, La, Ce, Pr, 
Nd, Sm). In order to study the s-process pattern in more detail, 
four indices are generally defined: [s/Fe], [ls/Fe], [hs/Fe] and 
[hs/ls]. Which elements are taken into account to determine 
these indices is different from author to author and is mainly 
determined by the possibility to compute accurate abundances 
of the different elements. To be able to compare our results with 
our results on the 21/j.m objects (Paper I), we define the ls-index 
as the mean of Y and Zr and the hs-index as the mean of Ba, 
La, Nd and Sm. Consequently, [s/Fe] is the mean of the former 
six elements and [hs/ls] = [hs/Fe] -[ls/Fe]. 

With the definitions from above, the s-process in- 
dices for IRAS 06530-0213 are: [s/Fe] = +2.1, [ls/Fe] = +1.8, 
[hs/Fe] = +2.2 and [hs/ls] = +0.4; for IRAS 08143-4406, 
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Fig.6. [el/Fe] values for IRAS 06530-0213 and IRAS 08143-4406. 
The errorbars are the line-to-line scatters; the dotted line represents 
[el/Fe] = 0. It is already clear from the figure that a higher neutron 
irradiation is expected for IRAS 06530-0213. 



they are: [s/Fe] = +1.5, [ls/Fe] = +1.5, [hs/Fe] = +1.5 and 
[hs/ls] = 0.0. Note that the Ba abundance for IRAS 06530- 
0213 had to be estimated in order to calculate the [hs/Fe] 
and [hs/ls] indices. We estimated this abundance in the same 
way as explained in Paper I, by using the tables of Malaney 
(1987) and an exponential distribution of neutron exposures 
of to = 0.4mbar _1 . The Ba abundance estimated from the La 
abundance is then [Ba/Fe] = +2.55. 

These abundance analyses make IRAS 06530-0213 and 
IRAS 08 143-4406 join the group of post-AGB stars that clearly 
display chemical evidence for a third dredge-up event. They 
share the same chemical signatures as the 21 /im objects dis- 
cussed in Paper I. IRAS 06530-0213 is s-process enhanced at 
approximately the same level as IRAS 05341+0852, the most 
s-process enriched intrinsic object found so far. 

4.6. Comparison with Hrivnak & Reddy (2003) 

An independent abundance analysis of IRAS 06530-0213 has 
recently been carried out by Hrivnak & Reddy (2003), based 
on high-resolution spectra taken with the 2.7 m telescope at 
McDonald Observatory. The authors derived a slightly lower 
temperature of 6900 K, resulting in a lower metallicity of 
[Fe/H] = -0.9. The abundances relative to iron of the elements 
that are in common in both analyses are, however, in good 
agreement. One should note that the analysis presented here 
is based on spectra with both a higher resolution (56 000 vs. 
45 000) and a higher signal-to-noise (100 vs. 30). As a conse- 
quence, our analysis does not only include more lines (149 vs. 
102), it also covers much more ions (28 vs. 17). 



Table 7. The metallicity and the s-process indices of the eight stars 
plotted in Fig. 7. Some values given in this table slightly differ from 
the values given in Paper I, due to minor revisions and updates that are 
described in detail in Reyniers (2002). 



object 



[Fe/H] [s/Fe] [ls/Fe] [hs/Fe] [hs/ls] 



IRAS 04296+3429 


-0.6 


1.5 


1.7 


1.5 


-0.2 


IRAS 05341+0852 


-0.7 


2.2 


2.0 


2.3 


0.3 


IRAS 06530-0213 


-0.5 


2.1 


1.8 


2.2 


0.4 


IRAS 07134+1005 


-1.0 


1.5 


1.6 


1.5 


-0.1 


IRAS 08143-4406 


-0.4 


1.5 


1.5 


1.5 


0.0 


IRAS 19500-1709 


-0.7 


1.2 


1.4 


1.0 


-0.4 


IRAS 22223+4327 


-0.3 


1.1 


1.4 


0.9 


-0.4 


IRAS 23304+6147 


-0.8 


1.6 


1.5 


1.6 


0.1 



5. Neutron exposure 

In Paper I (Sect. 6) we compared the strength of the neutron 
irradiation parameterized by the [hs/ls] index with both the to- 
tal s-process enrichment and the metallicity for the 21yum stars. 
A surprising result from this comparison was the strong corre- 
lation between the [hs/ls] index and the [s/Fe] index implying 
an increase in efficiency of the neutron nucleosynthesis with in- 
creasing third dredge-up efficiency. This relation was, however, 
hampered by the low number statistics on which it was based: 
only six stars. Here we have two more data points (Tab. 7), one 
of which having a very strong enrichment (IRAS 06530-0213). 

From Fig. 7 it is clear that the results of the two newly dis- 
covered objects strengthen the conclusions of Paper I. In the 
upper panel of this figure, the two data points for the newly 
discovered enriched stars nicely fit the correlation between the 
[hs/ls] index and the [s/Fe] index. The new correlation coeffi- 
cient for the eight stars is +0.96, the same as found in Paper I, 
while a least-squares fit gives [hs/ls] = 0.70[s/Fe] - 1.14. 
Especially the analysis of IRAS 06530-0213 yields an impor- 
tant data point to extend the relation towards stronger enrich- 
ments. In turn, the addition of the new objects dramatically il- 
lustrates the weak correlation between the metallicity and the 
[hs/ls] index {lower panel of Fig. 7). Whereas the correlation 
coefficient for these two quantities for the 21/im stars was still 
-0.55, it decreases to -0.19 after addition of the two new ob- 
jects. Mainly the strong enrichment of IRAS 06530-0213 com- 
pared to its only mild metal deficiency is a strong indication 
to suspect a large intrinsic spread of integrated neutron irradi- 
ations (or, alternatively, 13 C pocket efficiencies). One can even 
doubt the presence of any relation between [Fe/H] and [hs/ls] 
after inspection of Fig. 7. 

The correlation between [hs/ls] and [s/Fe] is obviously de- 
pendent on the specific choice of the elements involved in the 
calculation of both the [s/Fe] and [hs/ls] index. The choice to 
incorporate Y and Zr for the calculation of [ls/Fe] and Ba, La, 
Nd and Sm for [hs/Fe] was based on the choice made by Busso 
et al. (1995). Some of these elements (like Ba and Sm) are, 
however, not always observed and abundances of unobserved 
elements were estimated using the Malaney (1987) tables (see 
Sect. 4.5). It is an interesting exercise to investigate a possible 
dependence of the relation presented above upon the specific 
elements used in the indices. On Fig. 8, the relation between 
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Table 6. Abundance results for the two programme stars. Because of a time gap of 16 days between the DI2 and RED580 UVES setting for 
IRAS 08143-4406, we had to adopt 2 different sets of atmospheric parameters. The parameters for the RED580 setting are deduced using the 
iron lines. The parameters for the DI2 setting are deduced relative to the parameters of the RED580 setting, by a detailed study of 26 lines 
in the overlapping region. The following solar abundances are adopted: for the solar iron abundance we used the meteoritic iron abundance 
of 7.51; the references for the solar CNO are C: Biemont et al. 1993, N: Hibbert et al. 1991 and O: Biemont et al. 1991; for Mg and Si the 
latest Holweger values (Holweger 2001) were used; for La and Eu we took the recent values derived by Lawler et al. (2001a) and Lawler et al. 
(2001b) respectively; other solar abundances were taken from the review by Grevesse & Sauval (1998). Despite the fact that there are more 
recent values for some of the solar abundances (especially for the solar CNO), we take these references to ensure as much as possible the gf 
values that we have used in the present paper to be consistent with the adopted solar abundances. 
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[Fe/H] = -0.39 
DI2 setting 
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r eff = 7050K 
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ss: spectrum synthesis; bl: line is blended, but abundance is corrected for this blend, using a dedicated MOOG routine 



the neutron irradiation and the total s-process enrichment is re- 
visited with different elements for the indices. 

On panel (a) of this figure, the relation is drawn for the "de- 
fault" index definition. On panel (b), the Ce abundance is in- 



corporated in the 'hs' (and hence V) definition, but the points 
hardly change their position on the graph. On the middle panels 
the 'Is' index is defined by only one element, being Y on panel 
(c) and Zr on panel (d). This change has a rather large impact 
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Fig. 7. Comparison of IRAS 06530-0213 and IRAS 08143-4406 with 
the 21/im stars of Paper I. The values of the data points in this fig- 
ure are given in Tab. 7. upper panel: The correlation between the 
total enrichment of s-process elements ([s/Fe], which is the mean 
of [Y/Fe], [Zr/Fe], [Ba/Fe], [La/Fe], [Nd/Fe] and [Sm/Fe]) and the 
[hs/ls] index. The straight line gives the simple least-squares fit. lower 
panel: The [hs/ls] index as a function of the metallicity determined 
by the Fe abundance. The dotted line gives the simple least-squares fit 
[hs/ls] = -0.23 [Fe/H] -0.18. 

on the position of the points. The restriction to one abundance 
for the definition of the is' index does not seem favourable 
in this context. Especially if the is' index is solely defined 
by the Y abundance, the relation is clearly less well defined. 
In Sect. 6 we further focus on the problems with this particu- 
lar abundance. On the two lower panels, the 'hs' definition has 
been varied. On panel (e), the [hs/ls] definition is the same as in 
Busso et al. (2001). With this definition, only observed abun- 
dances are used, none were estimated. Moreover, the relation 
found in panel (a) is very well reproduced with this definition. 
Hence, restricting hs to La and Nd, the definition seems to be 
a very good alternative for the "default" definition. In the last 
panel, panel (f), a "minimal" definition is given with only two 
elements, Zr and La, which makes the scatter on the mean trend 
to become quite high, but still a strong hint to the proposed re- 
lation exists. 

6. Comparison with AGB nucleosynthetic models 

6.1. AGB stellar models 

Although there is general consensus that the s-process nucle- 
osynthesis is triggered by the engulfment of a small amount 
of protons of the hydrogen-rich envelope into the top layers 
of the He intershell (e.g. Busso et al. 1999, and references 
therein), there is not yet a selfconsistent nucleosynthetic AGB 
model. Different evolutionary codes use different assumptions, 
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Fig. 8. Study of the dependence of the proposed relation between the 
[s/Fe] index and [hs/ls] index upon the particular definition of the in- 
dices. Despite the scatter on the mean trend for some definitions, the 
relation holds for each definition presented here. Especially the defini- 
tion proposed by Busso et al. (2001) (panel e) seems to be an excellent 
alternative for the "default" definition of panel a. 

not only for the partial mixing zone, but also for the third 
dredge-up efficiency with or without overshoot assumptions 
(see comparative study of Lugaro et al. 2003, and references 
therein). Rotational mixing has recently been analysed and a 
strong effect seems to be a decrease in the neutron irradiation 
since the major neutron poison ( 14 N) is mixed into the 13 C-rich 
layers (Langer et al. 1999; Herwig et al. 2003). 

In this section we compare our observed s-process distribu- 
tions with the results of nucleosynthetic AGB models (Gallino 
et al. 1998; Busso et al. 1999, 2001). The nucleosynthesis pre- 
dictions are calculated by a postprocessing code grafted upon 
the FRANEC evolutionary models (Straniero et al. 1997, 2003) 
that span a grid of models with masses of 1.3 to 3.0 M and 
metallicities from solar to 1/20 solar. In the models, an ad-hoc 
13 C profile is introduced which follows basically an exponen- 
tial distribution (see Fig.l of Gallino et al. 1998) and which is 
burnt in radiative conditions during the interpulse phase, at a 
temperature of T ~ 0.9xl0 8 K. Although the physical origin of 
the 13 C pocket remains to be studied, the models are used here 
to estimate the strength of the 13 C pocket needed to explain 
the observed s-process abundances. The 13 C pocket choice la- 
beled "Standard Case" (ST) corresponds to 4xlO" 6 M of 13 C 
and it is named in this way because it has been shown to best 
reproduce the "main component" of the s-process in the solar 
system when applied to low mass AGB stars of [Fe/H] = -0.3 
(Arlandini et al. 1999). Other choices of 13 C pockets for indi- 
vidual AGB stars are then obtained by multiplication or divi- 
sion of this standard 13 C amount. As a consequence, different 
s-process patterns can be obtained for a fixed metallicity and 
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mass. A spread in [hs/ls] such as found in the lower panel of 
Fig. 7 is therefore naturally found in these models when vary- 
ing the strength of the 13 C pocket. 

The FRANEC models allow the third dredge-up mecha- 
nism to self-consistently occur after a limited number of pulses, 
using the Schwarzschild criterion and without invoking any 
extra-mixing. Other AGB models (e.g. Goriely & Mowlavi 
2000; Herwig et al. 2003) use different prescriptions. Note that 
in the FRANEC models the thermal pulses with third dredge- 
up stop prior to the end of the AGB evolution, when the en- 
velope is still relatively massive and the superwind has yet to 
develop. The mass loss prescription is also a critical point in 
each AGB model code, and different groups use different ap- 
proaches (see e.g. Blocker 1999, and references therein). The 
models presented here have been obtained using the mass loss 
parameterisation by Reimers (1975), with the choice rj = 0.3 for 
the free parameter. 

It is clear that there is not yet a full theoretical understand- 
ing of the AGB nucleosynthesis and third dredge-up phenom- 
ena, making the matching between the models and the observed 
s-process distributions to be interpreted with caution. We want 
to stress, however, that the s-process enriched post-AGB ob- 
jects are ideally suited to confront observed s-process distri- 
butions with theoretically predicted abundances, since (i) no 
dilution has to be accounted for (as in the case of extrinsic 
AGB-stars) and (ii) these objects fully completed their evolu- 
tion through the TP- AGB phase (contrary to genuine AGB stars 
from the M-MS-S-SC-C sequence). 

6.2. Model fits for the 21 stars of Paper I 

We made an attempt to fit the s-process distribution of the 21/j.m 
stars discussed in Paper I with the latest AGB models. The 
abundances were re-calculated using the latest MOOG version 
(April 2002). The abundances of IRAS 05341+0852 were re- 
determined using new VLT+UVES spectra. Although this new 
analysis is based on superior spectra (see Fig. 1) compared to 
the WHT+UES spectra in Paper I, the new results are not dras- 
tically different from our previous analysis. 

The model parameters were found in the following way. 
We restricted the models to an initial mass of 1.5 M . For the 
metallicity of the model, we used the spectroscopically derived 
value. The 13 C pocket is then the only parameter that is left 
over. It was found by manually fitting the observed Zr and 
La abundance (which is in fact representative for the [hs/ls] 
index). For two stars in the sample (IRAS 19500-1709 and 
IRAS 23304+6147) no solution could be found. For these two 
stars, the best fit was obtained with a model of lower initial 
mass of 1.3 M©. The fits are presented in Fig. 9. The fit for 
IRAS 19500-1709 still shows quite large discrepancies for the 
heavier elements. The abundances of these elements are, how- 
ever, derived from only a few lines (see Paper I, for details), 
and should therefore be considered with caution. 

The fits are not unique, but it turned out that a free choice 
between the possible 13 C pockets was rather well constrained 
due to the fixed metallicity. On the other hand, the amount 
of thermal pulses can influence the choice of the 13 C pocket: 
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Fig. 9. Abundances of the s-process elements for the six 21/im stars of 
Paper I, together with their best fit AGB model. The model parameters 
are given in the upper left corner; the initial mass is taken 1.5 M Q for 
all models, except for IRAS 19500-1709 and IRAS 23304+6147. For 
these two stars, a lower initial mass of 1.3 M Q yielded a better fit. The 
errorbars plotted on the figure are the line-to-line scatters. The actual 
error on an abundance is, however, in most cases considerably larger 
than this value, due to uncertainties in the atmospheric parameters, 
undetected blends, etc. Note that in the AGB model predictions the 
solar meteoritic abundances by Anders & Grevesse (1989) are used 
to convert the s-process yields to the actual [el/Fe] predictions. This 
specific choice has only a marginal effect on the predictions in this 
elemental range. Choosing the solar abundances as listed in Tab. 6 to 
calculate the predictions, would introduce a maximum change of only 
0.07 dex (Pr) in the predictions. 



e.g. Busso et al. (2001) propose a 13 C pocket "ST/3" with 10 
thermal pulses for the same abundances of IRAS 07134+1005 
we present here, whereas we found "ST/5". However, note that 
the present AGB model predictions have been considerably im- 
proved since the Busso et al. paper. 

Finally, a small remark is given on the niobium (Nb, Z = 41) 
abundance. Since these models are for intrinsic stars, 93 Zr with 
a half-life of 1 .5 x 1 6 yr, has not yet decayed to Nb. In extrinsic 
stars, the Nb abundance is expected to be at the same level as 
Zr. 
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Fig. 10. Abundances of IRAS 06530-0213 and IRAS 08143-4406 to- 
gether with their best fit AGB model. The model parameters are given 
in the upper left corner. The initial mass is 1.5 M Q . Errorbars are the 
line-to-line scatters that are given in Tab. 6. Note that the use of the 
solar meteoritic abundances in the AGB model predictions (see also 
the caption of Fig. 9), introduces a rather large uncertainty for the 
W abundance, for which the difference between solar photospheric 
and meteoritic abundance is very large (0.43 dex, Anders & Grevesse 
1989). 



6.3. Model fits for the two programme stars 

After this first exercise, we searched for the best model fit for 
the two programme stars of this paper. The same fitting method 
was applied as for the 21/im stars, except that we extended the 
comparison to all observed elements. The final model fits, to- 
gether with their parameters, are shown in Fig. 10. Note that 
decreasing the metallicity by e.g. 0.3 dex (which corresponds 
to the metallicity found for IRAS 06530-0213 by Hrivnak & 
Reddy 2003) would result in a fit with a 13 C pocket which is 
decreased with the same factor, with respect to the original one. 

6.4. Model fit conclusions 

Apart from the fact that several solutions are possible for the 
same pattern, it is clear, though, that only one 13 C pocket 
alone cannot explain the abundance patterns in all eight ob- 
jects. IRAS 07134+1005 for example can certainly not be fit- 
ted with a 13 C pocket which exceeds "ST/2", whereas for 
IRAS 05341+0852 a small pocket (<"ST") is not possible. This 
difference in 13 C pockets is much larger than the difference that 
could result from uncertainties in the models or in the abun- 
dances, and therefore we believe that the spread in 13 C pock- 
ets is real. The need of an intrinsic spread in the adopted 13 C 
pockets should be regarded as an essential result of this first 
exercise. 



Other interesting conclusions can be drawn from the fits by 
inspecting the abundances individually. We give four examples 
to illustrate the mutual benefit of such fits for both the models 
and the observational analysis. 

(1) The carbon abundance (not on the figure for the 2\ym 
stars) is predicted far too high in each fit. Since this particu- 
lar abundance is very reliable for each star (cr< 0.2), the too 
high prediction is probably a model's artefact. 

(2) The yttrium abundance (Z=39) is systematically 0.2 dex too 
high for the cooler objects (r e ff 7000 K) in the sample when 
compared to the predicted abundances. For these cooler stars, 
only a few lines (typically ~4) were found since most Y n lines 
are too strong to deduce an accurate abundance. The two hot- 
ter stars in the sample (IRAS 07134+1005 and IRAS 19500- 
1709) do not show this effect, and since the hfs effect for Y is 
only very small, this discrepancy between the observed and the 
predicted abundance for the cooler stars is probably due to a 
log(g/) problem and/or blending of the Yn-lines used in these 
stars. 

(3) Another, although minor, discrepancy is that the cerium 
abundance is systematically too low compared to the predic- 
tions for at least five objects. In this context it is noteworthy that 
the Ce abundance decreases by another ~0.2 dex when calcu- 
lated by using the log(g/) values of the D.R.E.A.M. database. 
Hence, replacing the VALD data by the D.R.E.A.M. data for 
Ce would only increase the observed discrepancy. 

(4) The predictions for the elements beyond the Ba-peak for 
IRAS 06530-0213 and IRAS 08143-4406 are very consistent 
with the observed abundances. This is an additional argu- 
ment strengthening the line identification for these elements. 
Moreover, their abundances are very sensitive to the adopted 
13 C pocket. Hence, these elements (especially Hf and W) are 
ideally suited to discriminate between possible 13 C pockets. 
Unfortunately, a Hf abundance is very difficult to derive in 
these objects since it has only suitable lines in the blue, and 
the W abundance derived from the 5 104.432 A line depends 
on the adopted Sm abundance (see Reyniers & Van Winckel 
2003). 



7. Conclusion 

IRAS 06530-0213 and IRAS 08143-4406 were classified as 
post-AGB candidates on the base of their IR properties and 
their spectral type. The analysis presented in this paper proved 
that both objects also show clear chemical evidence of their 
evolved, post third dredge-up character, since they show not 
only a clear carbon enhancement, but also a large enrichment in 
s-process elements. Especially IRAS 06530-0213 is an object 
that is very interesting for further research since it is enhanced 
at almost the same level as IRAS 05341+0852, the most s- 
process enriched object known so far. Moreover, detailed abun- 
dances of elements beyond the Ba peak were obtained, a result 
which was possible due to the combination of the high qual- 
ity VLT+UVES spectra and newly released line data in both 
the VALD and the D.R.E.A.M. databases. The end point of the 
s-process (Pb) is expected to be produced at important levels, 
also in these moderate metal-deficient objects, but Pb has no 
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suitable lines in the optical spectrum for these temperatures and 
gravities. 

The atmospheric and chemical properties of the two ob- 
jects are remarkably similar than those of the 21yum objects. 
Especially the relation between third dredge-up efficiency and 
neutron nucleosynthesis efficiency outlined in Paper I for the 
21yum stars is confirmed by these two new objects. On the other 
hand, the results on the two objects made the relation between 
metallicity and neutron nucleosynthesis efficiency even more 
confusing: whereas the 21yum stars of Paper I still displayed a 
weak correlation, the new data points (especially IRAS 06530- 
0213) suggest that there is no correlation at all. This means that 
there is a large intrinsic spread in neutron nucleosynthesis ef- 
ficiency. Such a spread points to different 13 C pocket strengths 
in probably highly similar objects which is confirmed by our 
modelling, in which quite different 13 C pockets are needed to 
model the s-process distribution. A physical explanation for the 
different 13 C pocket strengths is not yet found. The proton en- 
gulfment induces an interplay between the 14 N pocket, acting 
as a neutron poison, and the 13 C pocket, acting as a neutron 
donor. In this paper we restrengthen the observational finding 
that even for objects with similar metallicities, quite a different 
neutron irradiation is observed. Possibly rotation can play an 
important role in this discussion (Langer et al. 1999; Herwig 
et al. 2003). 

Moreover, apart from the s-process enriched objects dis- 
cussed in this paper, there are very similar post-AGB objects 
without showing any s-process enrichment. This dichotomy be- 
tween s-process enriched post-AGB stars, and post-AGB stars 
showing no s-enrichment at all, is also not understood. Note, 
however, that in the FRANEC models the third dredge-up does 
not occur below a certain critical mass, which slightly de- 
creases with decreasing metallicity. 

It is clear that a spectroscopic infrared study of the two ob- 
jects discussed in this paper would be an invaluable supple- 
ment to this optical study. Indeed, till now all enriched objects 
display the 21yum feature in their IR spectrum if this data is 
available. It should be tested if this assertion still holds after 
observing IRAS 06530-0213 and IRAS 08143-4406 in the IR 
spectral domain. 
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